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Abstract 

Given a static Schwarzschild spacetime of ADM mass Af, it is well- 
known that no ingoing causal geodesic starting in the outer domain 
r > 2M will cross the event horizon r = 2M in finite Schwarzschild 
time. In the present paper we show that in gravitational collapse of 
Vlasov matter this behaviour can be very different. We construct initial 
data for which a black hole forms and all matter crosses the event 
horizon as Schwarzschild time goes to infinity, and we show that this is 
a necessary condition for geodesic completeness of the event horizon. In 
addition to a careful analysis of the asymptotic behaviour of the matter 
characteristics our proof requires a new argument for global existence 
of solutions to the spherically symmetric Einstein- Vlasov system in an 
outer domain, since our initial data have non-compact support in the 
radial momentum variable and previous methods break down. 



* Support by the Institut Mittag-LefHer (Djursholm, Sweden) is gratefully acknowl- 
edged. 



1 



1 Introduction 



In a previous study [3j two classes of initial data for the spherically symmet- 
ric Einstein- Vlasov system were constructed which guarantee the formation 
of black holes. An additional argument to match the definition of a black 
hole in [8] is given [1]. In the present paper we denote any of these initial 
data classes by X. 

The analysis in [3] is carried out in Schwarzschild coordinates where the 
metric takes the form 



Here t G R is the time coordinate, r E [0, oo[ is the area radius, i.e., Anr'^ 
is the area of the orbit of the symmetry group S0(3) labeled by r, and the 
angles 9 £ [0, vr] and f G [0, 2tt] parameterize these orbits. The structure of 
the initial data X is such that a possibly large fraction of its ADM mass M 
is necessarily located in the outer domain r > 2M. In [3] it was shown that 
solutions launched by such initial data have the following property: there 
exist constants a, (3 > such that spacetime is vacuum for 



ds'^ = -[1 ] dt^+[l dr^ + r^{d9^ + sin^ 9d^^), 



representing a black hole of mass M. The generator of the event horizon 
approaches the surface r = 2M asymptotically as Schwarzschild time goes 
to infinity, cf. [3l Thm. 2.4]. 

Although (11.21) gives information about the asymptotic location of the 
matter it does not answer the question whether or not matter crosses the 
surface r = 2M. As a matter of fact, the inequality (II. 2|] is not sufficient 
to conclude that any matter initially in the region r > 2M ever crosses 
the surface r = 2M since matter can pile up at r = 2M. On the other 
hand it is known that not all matter can cross the surface r = 2M in finite 
Schwarzschild time. Indeed, if this were to happen the Einstein equations 
would imply that the metric function A became infinite at r = 2M. But 
according to [H] this cannot happen for the solutions considered in [3] . It fol- 
lows that on any finite time interval some matter must remain in the region 
r > 2M. The purpose of the present paper is to investigate the asymptotic 
behaviour of Vlasov matter in Schwarzschild time in the neighbourhood of 



ds^ 



(1.1) 



r>2M + ae"'^*, t > 0. 



(1.2) 



Hence in this domain the metric equals the Schwarzschild metric 
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the event horizon. Note that if matter crosses the surface r = 2M in finite 
time it also crosses the event horizon in finite time. Our main motivations 
are the following. 

• In Proposition 12.31 it is shown that a necessary condition for complete- 
ness of the outgoing radial null geodesic which generates the event 
horizon is that all matter crosses the surface r = 2M as Schwarzschild 
time goes to infinity. 

• In a static Schwarzschild spacetime of ADM mass M no ingoing causal 
geodesic starting in the outer domain r > 2M will cross the event 
horizon r = 2M in finite Schwarzschild time. It is interesting to know 
if this remains true in evolutionary gravitational collapse. The result 
in the present paper shows that for the initial data we construct the 
behaviour is indeed very different. 

• In [6l p. 13] some open problems about gravitational collapse are 
stated. For instance, for a scalar field it is known that 

supr = 2supm, (1-3) 

H H 

but for other matter models this issue is open. Here m is the quasi- 
local mass, H is the event horizon, and r the area radius. 

• [9l Thm. 1.5] relates the asymptotic behaviour of the matter at the 
event horizon to the question of strong cosmic censorship, see also [U 
Question 15.3]. 

• The asymptotic behaviour of matter in Schwarzschild time is directly 
related to what earth bound observers of gravitational collapse observe, 
which is not the case using other standard coordinates, e.g. Eddington- 
Finkelstein coordinates. If in Schwarzschild time all the matter crosses 
r = 2M, an earth bound observer will "see" all matter eventually 
swallowed by the emerging black hole. 

• An important open problem is the question whether or not solu- 
tions can break down in finite Schwarzschild time. It is often conjec- 
tured that Schwarzschild coordinates are singularity avoiding and that 
in these coordinates solutions of the spherically symmetric Einstein- 
Vlasov system exist globally for general initial data. We expect that to 
understand the asymptotics of Vlasov matter in gravitational collapse 
is going to be useful for understanding the global existence issue in 
general. 
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In the present paper we construct a class J of initial data for which the 
results in [3] apply, and such that all the matter asymptotically crosses the 
surface r = 2M. In particular (jl.3p holds for Vlasov matter for this class 
of initial data. The class J is different from the class X in that the support 
of the momentum variables is not compact. This is a technical condition 
needed for our method of proof, but we believe that the conclusion holds 
for compactly supported data as well. Below, we always have in mind the 
radial momentum variable when we discuss compactly or non-compactly 
supported initial data. Our method of proof does imply that matter crosses 
the surface r = 2M also in the compactly supported case, but we are not able 
to conclude that all matter eventually crosses r = 2M in this case. We point 
out that the condition of non-compact support is required in some works in 
the cosmological case, cf. [8] and [16]. For compactly supported initial data 
the result in [T3] guarantees that solutions exist as long as matter stays in a 
region r > e > 0. The proof in [Hj breaks down for non-compactly supported 
data. But for applying the method in [3j it is crucial that solutions are global 
in an outer domain, and so we need to establish such a global existence 
result in the case of non-compactly supported data as well. Non-compactly 
supported data have been considered for other kinetic equations such as 
the Vlasov- Poisson system [IT] and the Vlasov-Maxwell system |10^ I15j . 
These methods do not directly apply in the case of the Einstein- Vlasov 
system, and we have not been able to find a result analogous to [H] for 
non-compactly supported data. However, the initial data set J constructed 
below is such that the matter continues to move inward for all times. This 
crucial feature allows us to obtain the necessary global existence proof for 
the corresponding, non-compactly supported data in an outer domain. 

The outline of the paper is as follows. In the next section we introduce 
the Einstein- Vlasov system, recall the set up and the construction of the 
class of initial data in [3], and formulate the main results of the present 
paper. Sections 4, 5, and 6 are devoted to their proofs. 

2 Set up and main results 

In this section we recall the Einstein- Vlasov system and the set up in [3] 
and formulate the main results. For more background on kinetic theory and 
the Einstein- Vlasov system we refer to [1]. We consider the asymptotically 
flat spherically symmetric Einstein- Vlasov system. We use Schwarzschild 
coordinates {t,r,6,ip) and parameterize the metric as in (II. ip . Asymptotic 
flatness means that the metric quantities A and fi have to satisfy the bound- 
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ary conditions 

lim A(t,r) = lim /i(t,r) = 0. (2.1) 

r— ►oo r— >oo 

Vlasov matter is a collisionless ensemble of particles which is described by 
a density function / on phase space. In order to exploit the symmetry it is 
useful to introduce non-canonical variables on momentum space and write 
/ = f{t, r, w, L). The variables u; g] — oo, oo[ and L £ [0, oo[ can be thought 
of as the radial component of the momentum and the square of the angular 
momentum respectively. 

The Vlasov equation is given by 

dtf + e^-^|5,/ - (^Xtw + e^-Vr^ - e''"^;^) = (2-2) 
where 



E = E{r, w, L) := vT+^w^ipL/VS, 

and where subscripts indicate partial derivatives. The Einstein equations 
read 

e"2^(2rAr - 1) + 1 = SvrrV, (2.3) 
e"^^(2r/i,. + 1) - 1 = Svrr^p, (2.4) 
At = -Anrj, (2.5) 
and the matter quantities are given by 

^ poo /"OO 

p{t,r) = —l I Ef{t,r,w,L)dLdw, (2.6) 

^ J-OD Jo 

P(^^'') = ^ / -^fit,r,w,L)dLdw, (2.7) 
j{t,r) = - / wfit,r,w,L)dLdw. (2.8) 

^ J-oo Jo 

The equations (|2.2p - (|2.8p constitute the spherically symmetric Einstein- 
Vlasov system in Schwarzschild coordinates. For a detailed derivation of 
this system we refer to p2]- 

o 

As initial data we need to prescribe an initial distribution function / = 
/(r, w,L) > such that 

rr pr poo poo 

I AT:rfp{ri)dri = A-K^ I / / Ef{r],w,L)dLdwdri<-. (2.9) 
Jo Jo J-oo Jo 2 
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Here we denote by p the energy density induced by the initial distribution 
function /. If in addition the initial data is we say that it is regular. 
In previous investigations the condition of compact support was included in 
the definition of regular data, but compact support in w is not required in 
the present paper and is replaced by a suitable fall-off condition, cf. (j3.3|) 
below. The Cauchy problem is well defined for regular initial data. We will 
restrict ourselves to a smaller class of regular initial data which guarantee 
the formation of black holes. Clearly, black holes do not form for any initial 
data, e.g., if the data is sufficiently small matter disperses and spacetime is 
geodesically complete, cf. [13]. 

Let us recall the set up and the properties of one of the initial data sets 
constructed in [3j. We fix < ro < ri, and let 7+ be the outgoing radial 
null geodesic originating from r = rg, i.e., 



ds 



.(s) = e{/^-A){.,7+(«))^ ^+(0) = ro. (2.10) 



We consider solutions of the spherically symmetric Einstein- Vlasov system 
(|2.2p - (|2.8p on the outer region 

D := {{t, r) G [0, oo[V > 7^ (*)}• (2-11) 

Note that characteristics of the Vlasov equation can pass from the region 
D into the region {r < 7^(t)} but not the other way around so that initial 

o 

data / posed for r > tq completely determine the solution on D. 

Let M := ri/2 be the total ADM mass and define the quasi- local mass 

by 

;>oo 

m(t,r) = M-ATT p{t,r])rfdr]. (2.12) 

J r 

Let Mout < M be given and such that 



2(M - M< 



out J 



Take Ri > ri such that 



and define 



n - ro 



< n- (2-13) 



Ri-ri< 



Ro := ^in + Ri). 



We require that all the matter in the region [ro, oo[ is initially located in the 
strip [Rq, Ri], with Mout being the corresponding fraction of the ADM mass 
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M, i.e., 

foo pRi 

47rr^/5(r) dr = Aur"^ p{r) dr = Mout- 
/ro J Ro 

Furthermore, the remaining fraction M — Mout should be initially located 
within the ball of area radius ro, i.e., 



pro 

/ Anr'^pir) dr = M - Mo 
Jo 



'0 

If one considers the Einstein- Vlasov system on the whole spacetime r > 0, 
then the definition (j2.12p for the quasi-local mass is equivalent to the more 
standard one, namely m{t, r) = Air p{t, rj) rfdr]. This is because the ADM 
mass M = m(t, oo) is conserved. On the outer domain D the definition 
(j2.12p is more suitable, since it does not refer to the matter inside {r < 
7~''(t)} except for the fact that this matter is there and contributes to the 
total mass. 

The properties above concern the structure of the initial data in space. 
We also need to specify conditions on the momentum variables. Let W- < 
and Li > be given. In f3j we introduced the general support condition: 

o 

For all (r, w,L) G supp /, 



and if r G [Rq,Ri] then 



w < 0< L< L 



and 



3L 

< L < — rf)(r]) + T] rh{ri)^ r] £ [vq, Ri]. 
V 



Here we use the notation rh when p = p in (|2.12p . In addition to the 
conditions above the initial data Z in [3] were assumed to have compact 
support. If W- is sufficiently negative, a black hole of ADM mass M forms 
and limj^oo 7* (t) = 2M for a certain radially outgoing null geodesic which 
is the generator of the event horizon, cf. [3, Thm. 2.4] and [H Sect. 4.3]. The 
initial data we construct below do have the properties specified above, but 
the support in the radial momentum variable w is not compact. However, 
the proof in [3j goes through unchanged also for such data provided the 
solutions are global on the domain D. With respect to global existence in 
the domain D the following holds. 
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Theorem 2.1 Let regular initial data f be given with the properties speci- 
fied above, and such that the fall- off condition \3. 3|) is satisfied. Then the 
corresponding solutions of the spherically symmetric Einstein- Vlasov system 
Ii2. g)) -< [Oj) in the domain D exist for all t >0. 

We can now state the main result of the present paper. 

Theorem 2.2 There exists a class of regular initial data for the spheri- 
cally symmetric Einstein- Vlasov system such that the corresponding solu- 
tions have the asymptotic property that 

hm m(t,2M) = hm m{t,j*{t)) = M. (2.14) 

t—^00 t^OO 

As mentioned in the introduction the condition ()2.14p is a necessary condi- 
tion for completeness of the generator 7* of the event horizon. We state this 
in a proposition. 

Proposition 2.3 A necessary condition for future completeness of the gen- 
erator 7* of the event horizon is that \2.14^ holds. 



3 Proof of Theorem EOl 

Compactly supported, regular initial data launch a local regular solution 
which can be extended as long as the momentum support of the solution 
can be controlled \i2\ I13j . We do not give a complete proof for the cor- 
responding result for non-compactly supported data and restrict ourselves 
to establishing the main a-priori bounds. To this end it is convenient to 
introduce the Cartesian coordinates x = r(sin0 cos 99, sin sin cos 0) G 
with corresponding momentum variable S so that 

X ' V ^ I i9ii9 9 E /.-.-,\ 

w = , L = \x X vf, \vf = + —. (3.1) 

Here • denotes the Euclidean scalar product and \v\ the induced norm. We 
denote by {X,V){s,t, x,v) the solution of the characteristic system of the 
Vlasov equation, written in the variables x and v, 

v = - (^At(s, r) ^ + e(^-^)(^'^)/i^(s, r) y^T+W) 
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with {X,V){t,t,x,v) = {x,v); here ' denotes the derivative with respect to 
s. We define 

Q{t) := sup I ,^1^(0^';^'^^^), \0<s<t, (X, .) e supp /(.) 

= sup| ^ + '^(_^'°;^"'"^' \0<s<t, (X,.) G supp/} . (3.2) 

We require that the initial data satisfy the fall-off condition 

11/11 := sup {l + \v\f\fix,v)\<oo. (3.3) 

f{t,x,v) = f{{X,V){0,t,x,v)), 

we get the estimate 

f{t,x,v) < 11/11(1 + \v{o,t,x,v)\r' < QHrnim + \v\r'. 

Hence 

{l + \v\)f{t,x,v)dv<\\f\\Q^{t) [ il + \v\)-Uv<C\\f\\Q'it). (3.4) 



By the characteristic system, 

^(1 + \V{s,0,x,v)\) < (||e^-V.(^)l|oo + ||At(s)|U)(l + \V{s)\). 
This implies that for < s < t, 

l + \Vis,0,X,v)\ ^ g/„*(||e''-V.W||.o + ||A,(r)||^)dr 
l + l^l - 

and we obtain the estimate 



The field equations ()2.3p and (j2.4p together with the boundary condition 
(I2TD imply that 

POO 

(/i + A)(t,r) = - / {fir + K){t,r])dr] <0, 
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and 



e'^-Vr(t,r) = e^+^ ( + 47rrp{t , r) ] < 47rr (| U + ||p(t)||c 



Together with (12. 5|) and (13. 4|) we have 

||e^-Vr(s)||oo + ||At(s)||oo < C(l + s)||/||Q'(s), 

so that 

This impHes that Q is bounded on some time interval [0,T[. A standard 
iterative procedure then shows that there is a local, regular solution which 
can be extended as long as the function Q does not blow up, cf. [T2l [T3] . 

Global existence in the outer domain D will now follow if we can establish 
a bound on Q{t) in D. For this argument we use the variables (r, w, L). By 

(1331), 

f{r,w,L) < C\w\-^. 



Let us define a quantity as in ()3.2p . By abuse of notation we let 
f 1^1 

^^'^'■= 1 \W{0,s,r,w,L)\ I ^ ^ ^ ^' ^) ^ ^^PP 

= sup|&l^jll^|0<.<t, (r,z.,L)Gsupp/}; (3.5) 

notice that in the outer domain D the area radius r > rg > so that by 
()3.ip a bound on Q as defined in (|3.5p implies a bound on Q as defined 
in (|3.2p . The following lemma taken from [3j shows that when the general 
support condition holds, then the particles in the outer domain D keep 
moving inward in a controlled way. 

o 

Lemma 3.1 Let f be regular and satisfy the general support condition for 
some suitable < 0. Then for all characteristics {R{t),W{t), L) with 
(i?(0), WiO), L) G supp / and R{0) € [Rq, Ri], 

W{t) < eA(t,iJW)g-A(o,i?(o))^(o) < e-A(o,i?(o))^(0) 

as long as (t,R{t)) £ D. In particular w < for all {r,w,L) G supp/(t) 
and {t, r) G D, and j <0 on D. 
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Following [13] we find that along any characteristic in supp/, 

J r-QO i-Li 

—w^<Cw'^ + Cj I \w\f{s,r,w,L)dLdw; (3.6) 



oo JO 

for this estimate it is essential that all particles are moving inward. We 
estimate the last term. Since 

o 

f{s, r, w, L) = f{RiO, s, r, w, L), W{0, s, r, w, L), L), 

we have 

oo /-Li 

/ \w\f{s,r,w,L)dLdw 
-oo Jo 

<C ff \w\\W{0,s,r,w,L)\-^dLdw. 

J Jsupp/(s,r,-,-) 

By the definition of Q and the general support condition, 

\w\ 



\WiO,s,r,w,L)\ >max<^ 
Hence we find that 

f°° /"^i ~ ~ /.oo |.Li / J- w ]\'^ 

/ / \w\f(s,r,w,L) dLdw < C / w { max < \ W^\, > dw 
J-ooJo Jo Jo V I Q[s)}J 



\W-\Qis) r-oo Q3fg\ 

< C I w dw + w :T—dw 

lo J\W.\Qis) W-^ 



< CQ\s). 
We have thus derived the estimate 

< Cw^ + CQ^(s), 

ds 



and hence for < s < t, 

.,2 



+ C J^Q\T)dT 



This implies that 



Q^(t) <e^* (^1 + C Q\s) ds^ 



hence Q is bounded on bounded time intervals, and global existence in D 
follows. □ 
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4 Proof of Theorem D 

We aim to show that ah characteristics starting in the domain [Rq , Ri [ enter 
the region {r < 2M} in finite time, and we need an estimate for the required 
time. Let {R{s),W{s), L) be a characteristic emanating from the support 

o 

of / with i?(0) G [Rq,Ri]; all the fohowing estimates are vahd as long as 
R{s) > 2M. By Lemma |3. II and the characteristic equation, 

where 

nir,w) : = 



^l + e-2M0,r)y;2 + _^^(2M)-2 

o 

We require that on supp /, 

uj < -e^(°''^)i^(r) (4.1) 

where K : [i?o, -Ri[^]0, oo[ is an increasing function which wih be specified 
befow. Hence 

K(r) 

B(r,w) > ^ ' := B(r), 

- y/l + K^{r) + Li{2M)'^ 

and 

R{s) < -B{R{0)) e(A'-^)(«-«{^)) . (4.2) 
By [21 Lemma 4.1 (a)], 

/Li — A > 2/i 

where 

Jr T 

Inserting this into ()4.2p imphes that 

R{s) < -B{R{0)) e^/^^^'^W) . (4.3) 

In order to estimate the right hand side we compute fit, cf. [31 Lemma 4.1 
(d)], and observe that 



fit{s,r)= / 47rj(s,??)e(^+^)(^''')e2^("''')(ir/ 

Jr 
1 f°° 

>— 47rj? 2j{s, rj) e('^+^)(^' ^^e^^^^' ^) dr/; 

2^" Jr 
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note that by Lemma [3?H j < 0. Since 



2 / \ 2 

E + 2w + — = { VE + ^] >0 



the expressions for the matter terms imply that 2j > —{p + p) so that by 
P Lemma 4.2], 



1 f °° 

Moreover, 



- . . m(s,r) 
p-r{s,r) = 2 — ^ 



,2A(s,r) 



r 



and |^(s)| < e(/^-^)(^'^(«)). Hence 



-/i(0,i?(0)) = / — fl(s,R(s))ds 
Jo as 

fit{s, R{s)) + fir{s, R{s))R{s)) ds 



- ' ' 2R{s) V ) R{sf ^ ) 

The right hand side of this inequahty will be estimated using the following 
lemma. 

Lemma 4.1 Let p € L^([2M, i?i]) he such that p>0 andO < 2m{r)/r < 1, 
where 

rRi 

m{r) := M — A.-nrf p{'q) dr], 

J r 

and let 

e^AW ,= rG[2M,i?i]. 
Then for all r G [2M,Ri], 

^-2/«i4^r,p{,,)e2M^)d,, - r - 2M 



> 



r — 2m{r) 
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Proof. For r G [2M, Ri] we define 

/i(r) := (r - 2m(r)) e^^'') - r + 2M, ^(r) := -2 / ATrr]p{r])e'^^'-''Ur]. 



Then 

/i'(r) = (1 - 8TTr^p{r)) e^W + (r - 2m(r)) e9W5'(r) - 1 

= (1 - 87rrV(r)) e^^"^) + r (^1 - '^H^^ e^'('')8^rp(r)e2^(^) - 1 

= (1 - 8irr^p{r)) e5(^) + re-2^(^')e^'('')87rr/)(r) e^^^''^ - 1 
= e^^*^) - 1 < 0, re [2M,i?i]. 

Hence for r G [2M,iii], 

/i(r) > h{Ri) = 0, 

which is the assertion. □ 

Remark. It is interesting to note that the configuration for which equality 
holds in the inequality in the lemma can be shown to be an infinitely thin 
shell. This should be compared to the situation considered in j2j where an 
infinitely thin shell is the maximizer of a similar integral expression as above. 

Let us return to the proof of Theorem 12.21 By the field equations (j2.3p and 
(12. 4p and the form ()2.6p and ()2.7p of the matter terms, 

+ Ar = Airre^^ip + p) < Bvrre^^p. 

Hence, Lemma l4.ll implies that for r £ [2M, 

(^+X){s,r) > p-2//'l 47rr?p(s,,y)e2M=,^)dr, - r - 2M 



r — 2m{s, r) 

We insert this into the estimate (j4.4p and find that as long as R{t) > 2M, 

p{t,R{t)) > /i(0,i?(0)) - J^ds. 
By ()4.3p this implies that 

R{s) < -B{R{0))C{Rmexp (^_2£-^cir) , 
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where C(r) := e'^'^^'^'^\ This imphes that 
d 1 _ R{s) 

which upon integration yields the estimate 

1 1 BC f ( M / 



R{t) - R{0) 2M V V Jo 

- RiO) ^ 2M V )■ ^ 

In order to proceed the functions C and B must be related properly. We 
require that 

p{r) > for r e]Ro,Ri[ 
so that rn(r) < M for r G [0, Ri [, and 

/ 2M , \ 2M 

We can therefore choose the function K which specifies our support condition 
(|4.1|) in such a way that for r G [i?0)-Ri[) 

K(r) 2M 

B{r) C{r) = ^ ' C{r) > 1 ; (4.6) 

^ ^ ^ ^ ^l + K2(r) + Li(2M)-2 ^ ^ r ' ^ ^ 

note that this necessarily implies that K(r) ^ 00 as r — > and that the 
function C is determined by the initial data. Given any r* £]Rq,Ri[ there 
now exists k> such that for r £ [RQ,r*], 

B{r) C{r) >1- — + K. 

r 

The estimate ()4.5p therefore implies that for any characteristic as above, but 
with i?(0) G [i?o,r*], and as long as R{t) > 2M, 

1,1 1 ^. 2M \ -2Mt/Rl\ 



- R{0) 2M V i?(0) 



2M V y i2(0) 

> ^ + ^ /"i _ ^-2Mt/Rf\ 
2M \ )' 
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which imphes that 



R{t) <^(l- e-2^^*/«? 

^ ^ 1 + K V 



This shows that there is a time t* > such that R{t*) < 2M for all charac- 
teristics starting with -R(O) G [Ro,r*]. 

To complete the proof we fix e > and let < Ri be sufficiently close 
to Ri such that rn(re) > M — e. Then there is a finite time such that all 
characteristics {R{t),W{t), L) with R{0) G [RQ,re] reach r = 2M at some 
time t < t^. 

We construct a curve (t, a(t)) with the property that m{t, a{t)) > M — 2e 
for <t <te and a{t) = 2M for some t < t^. To this end, let 

(4.7) 



ATTRjt, 

and let a be the solution of 

«= ''l''"i7% (^-^)(-'"\ a(0) = r,. (4.8) 
/3(s, a) + d 

The reason for introducing the 6 parameter is to avoid any potential prob- 
lems with uniqueness of solutions if p = 0. Taking the partial deriva- 
tive of e~^'^ = 1 — 2m /r with respect to t and using ()2.5p we find that 
rrit = —AiTr'^e^~'^j, and hence 

^m{t,a{t)) = -ATre'^^'-^'>'^''"^''>^a'^{t)j{t,a{t))+4TTa\t)p{t,a{t))a{^ 



= 4TO^e^'^-'^^^''°WJ ± 

> -47ra2(t)e(^-^)(*'"W)(5> --. 

Here we used that j < 0. Thus for all < t < t^, 

m{t,a{t)) > M - 2e, 

and it remains to show that a{t) < 2M for some < t < ■ Define for each 
r e [Ro,Ri] the barrier curve {t,RB{t)) by 

Rb = -5(r)e(^-^)(^'^«), iiB(O) = r. 

We use the term barrier curve since the area radius along this curve is larger 
than the area radius along any characteristic {R{t), W{t),L) starting in the 
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support of / with R{0) = r. This is clear from the differential estimate (j4.2p 
and the estimates which followed. In addition, all barrier curves starting at 
some r E [RQ,rf] reach the region r < 2M within the time interval [0,te]. 
For r G [Ro,Ri] the definition of B{r) and the condition on the support of 

o 

/ imply that, 

poo poo I 

|i(0,r)| > e^(°'")s(r)^ / / Jl + e-^^i^'^^w^ + Li{2M)-^ fdL dw 

f J-oo Jo 

>B{r)p{0,r), 



and since S < 1 we thus have 

\j{0,r)\+6 



> B{r). 



p(0,r)+5 

Consider the barrier curve Rsit) with Rb{0) = r^. Then 
p(0,re) +5 

and hence a{s) < -R_b(s) on a time interval ]0,,si]. Assume that a{t) > 2M 
for all < t < t,. We define 

t* := inf {t G [0,te] | -R_B < « on [t,te] for all barriers starting in [i?o,re]} . 

Then there exists some barrier curve Rsit) starting at some r G [iioj'^e] 
such that a(t*) = r* := Rsit*). By definition of t*, a{t) > Rsit) for t > t*, 
and hence 



p{t*,r*) + J 
which implies that 

The latter inequality is only possible if there is at least one characteristic 
{R{t),W{t),L) with R{t*) = r*, and 



Vl + i^2(^^(0)) + Li(2M)-2 ^l + |H^(t*)|2 + Li?(t*)-2 

X(/2(0)) 



> S(i?(0),W^(0)) > S(i?(0)) 



Vl + ^2(^(0)) + Li(2M)-2' 
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Since the function K is taken to be increasing this estimate imphes that 
R{0) < Rb{0) < r,. The barrier curve {t,RB{t)) with ^^(0) = R{0), must 
satisfy the estimate a{t*) = r* = R{t*) < Rsit*), and this is a contradiction 
to the definition of t*. Hence a{t) < 2M for some < t < t^, which proves 
that hmt^oo m{t, 2M) = M. 

We can chose the parameter 5 in (j4.7p such that m{t,a{t)) > M — 2e 
for all < t < + 1. Since a is strictly decreasing, a{te + 1) < 2M, 
and since limt^oo7*(i) = 2M there exists some time t > + 1 such that 
l*{t) > o(ie + 1) and hence m(t, 7*(i)) > M — 2e. This completes the proof 
of Theorem [221 □ 

To conclude our main result we show that initial data which satisfy the 
conditions required above do exist. To this end, let p = p{r) be a function 
supported in [Rq,Ri] with p{r) > on ]Rq,Ri[, and such that 

Mont = 47r / r'^p{r)dr < M, 
Jro 

satisfies (I2.13p . Define 

m(r) := M - 4^ H r^p{r)dr, e'^^^^^ := 1 - 



If r G [Rq,Ri[, then M — m{r) > 0, and hence there exists a function K 
as introduced in (|4.1|) . which satisfies the condition (j4.6|) . Now let h = 
h{r,w,L) be a function supported in [i?Oi-Ri]x] — oo, oo[x [0, Li] and 
such that h{r,w,L) = if e~^^'^^w > —K{r), and 



h{r, w, L) dL dw 



Let 

h{r, w, L) 



h{r, w, L) 



so that p{r)h{r, w, L) induces the energy density p and the quasi-local mass 
m. Let fi be a density function supported in [0, rg] such that the assumptions 

o 

of Lemma [3 . 1 1 hold . Then /= fi+ ph defines initial data which have all the 
required properties stated above. 

5 Proof of Proposition 12.31 

The function m(t, r) is increasing in both variables. We assume that 
m(i,r) < m(t,2M) < C < M for all t > and r < 2M and have to 
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show that the generator 7* of the event horizon is incomplete. According to 
[21 Thm. 2.4], 7* approaches r = 2M as t ^ 00. It follows that for ah t > 0, 

X{t,j*it))<C. (5.1) 

Let r I— > {t, r, 9, (f)){T) be an affine parameterization of 7* with corresponding 
momenta {p^ ,p'^ ,p^){t). Since 7* is radial, let 6 = Tr/2, = and 
p2 _ p3 — Q_ Si]2ce 7* is null, 

and we get 

pO = e^-V- (5.2) 

By the geodesic equations, 

^ = -e'(^-'^ prip'? - Kip'f - 2X,pV 
= 47r7*e2\pi)2[-p-p + 2j]. 

Here we used (|5.2|) to express p^ in terms of p^ . Since dt /dr = p^ , 

4^7*e'^+V[-(p + /5) + 2j], (5.3) 



dp^ 



and 7* is incomplete if 



By ^ and 



""^ ds 

< 00. 



r°° ds 

Jo P° 



'''' \s = /* ef^-^^e^o^^re^+Hp+P-Wdvas. (5.4) 



/o P° 7o Jo P^O) 

For the interior integral expression we have since p < p, 

/ 4^7*e^+^(p + p-2j)dr < / 4^^-f*e^'+^{2p - 2j)dT. 
Jo Jo 

Along the null geodesic we get 

^A(.,7*(.)) = A, + A.^ = 4.7* (p - j)e^^' - j^e^^\ 
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Hence by (f5l^ . 

47r7*e'^+^(2p-2i)dr = 2A(s, 7*(s)) - 2A(0, 7*(0)) 

^ r 2m(r,r(T)) ^(,,+^)(,^^.(,))^^ 

<C + C [ e^(^'^*(^))(ir. 

JO 

By [21 Thm. 2.4] there exist positive constants a and /3 such that if {t, r) 
satisfies 

r > 2M + ae"'^* := a{t), 

then there is vacumm at (t, r). By monotonicity of /z and the fact that there 
is vacumm for r > a, 

and thus ^ 

/ 47r7*e'^+^(2p-2j)(ir < C7. 

JO 

From (j5.4p we therefore obtain the estimate 

f%<cf e('^-^)(^'^*(^))d. = C f ^^ds < 2CM, 
Jo Jo Jo ds 

which says that 7* is incomplete, and the proof is complete. □ 
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